Phase diffusion and fractional Shapiro steps in superconducting quantum point 

contacts 
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We study the influence of classical phase diffusion on the fractional Shapiro steps in resistively 
shunted superconducting quantum point contacts. The problem is mapped onto a Smoluchowski 
equation with a time dependent potential. A numerical solution for the probability density of 
the phase difference between the leads gives access to the mean current and the mean voltage 
across the contact. Analytical solutions are derived in some limiting cases. We find that the 
effect of temperature is stronger on fractional than on integer steps, in accordance with preliminary 
experimental findings. We further extend the analysis to a more general environment including two 
resistances and a finite capacitance. 
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I. INTRODUCTION 

Shapiro stepsi in superconducting tunnel junctions are 
a clear evidence of the a. c. Josephson effect^, i.e. the flow 
of an alternating current through the junction under a 
finite dc voltage. Shapiro steps appear as a consequence 
of the beating between an applied microwave field and 
this alternating current. Their observation thus provides 
an indirect test of the sinusoidal current-phase relation. 
For the case of highly transmissive junctions, like atomic 
size contacts, the presence of higher harmonics in the 
current-phase relation leads to fractional Shapiro steps^. 

In contrast to conventional tunnel junctions, atomic 
contacts are characterized by a reduced set of conduc- 
tion channels whose transmissions r n can take arbitrary 
values between and I—. The channel content of a given 
contact can be determined experimentally with high ac- 
curacy by analyzing the subgap structure of the TV char- 
acteristic in the superconducting stated. Moreover, the 
electromagnetic environment of the contact can be de- 
signed by means of litographic techniques^. Due to all 
these properties atomic contacts can be considered as 
ideal systems to test theoretical predictions on meso- 
scopic electron transport under controlled conditions^. 

Shapiro steps in a superconducting atomic contact of 
arbitrary transmission were analyzed in Ref. within a 
fully microscopic approach in which an ideal voltage bias 
on the contact was assumed. However, a more realistic 
description of this phenomenon requires to take into ac- 
count environmental effects, the most important being 
certainly the effect of phase diffusion caused by thermal 
noise in the circuit in which the contact is embedded. 

The influence of thermal fluctuations on the current- 
voltage characteristics in tunnel junctions has been tradi- 
tionally analyzed by means of the so-called resistively and 
capacitivcly shunted junction (RCSJ) models. The start- 
ing point of this approach is to write down a Langevin 
equation for the phase difference accross the junction^. 
In this article, we generalize this model to the case of 
highly transmissive junctions in the presence of a mi- 



crowave field in order to study the effect of temperature 
on fractional Shapiro steps in atomic size contacts. Our 
work is also motivated by experiments underway in the 
Quantronics group at the C.E.A. in Saclay. Their pre- 
liminary results seem to indicate that fractional steps are 
more affected by thermal fluctuations than the integer 
onesiS. 

This study will be based on two main assumptions: 1) 
the external frequencies and thus the range of voltages 
considered are small compared to the superconducting 
gap on the contact leads. This would allow us to assume 
an adiabatic response of the contact following its static 
current-phase relation and 2) the resistances in the circuit 
containing the contact are small compared to the resis- 
tance quantum so as to neglect quantum effects in our 
treatment. As we discuss below, low values of the shunt- 
ing resistances are also necessary in order to observe well 
defined Shapiro steps even when thermal fluctuations are 
negligible. 

The paper is organized as follows: In the first sec- 
tion we present our generalization of the RCSJ model 
for contacts of arbitrary transmission in the presence of 
microwave radiation. The second section describes how 
to solve the stochastic equation of the problem and shows 
the numerical results for the current-voltage characteris- 
tics. In the third section we discuss some limiting cases in 
which analytical results can be obtained. In the last sec- 
tion we consider a more general electromagnetic environ- 
ment including two resistances and a finite capacitance 
and discuss its effect on the supercurrent peak in the ab- 
sence of radiation. We finally present some concluding 
remarks. 



II. THE MODEL 

A. Electrical circuit 

We consider the standard RCSJ model, to which we 
add a microwave field. The equivalent circuit is shown on 
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FIG. 1: Equivalent circuit of the resistively and capacitively 
shunted junction (RCSJ) model in presence of a microwave 

field Vac COS Ldt. 



the figure^] The parallel combination of the contact with 
a resistance R and a capacitance C is current polarized 
by a current source The microwave field introduces 
an additional ac voltage, v ac on the branch containing 
the contact. We denote by 9 the superconducting phase 
difference across the contact and by 1(9) the correspond- 
ing current-phase relation. Dissipation in the resistance 
will cause Johnson- Nyquist noise L(t). Conservation of 
the current implies that 



(1) 



with V = v + v ac cos uit. The voltage across the contact 
is related to the phase via the usual relation v — 4>q9 
(0o = h/2e being the reduced flux quantum). 

In order to stress the analogy with the Brownian mo- 
tion of a particle in a time-dependent potential, we fol- 
low Ref. M and define M = $C, p = <j) vC, 77 = 1/RC, 
r(t) ee -<f> L(t)/r]M and 
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The originality of this work lies in the explicit time de- 
pendence of the potential U(9, t) and the non-trivial form 
of the current-phase relation 1(9) which we describe in 
more detail below. When the capacitance is negligible, 
we are in the strong damping regime r\p ^> p, and JTJ 
takes the form of a simple Langevin equation 
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We consider a gaussian white noise, that is 
< r(t) > = 



< r(t)r(t') > = S(t' - t) 
V[r(t)} 



exp 



2T 
77M 

~4T 



r 2 (t)dt 



(4) 
(5) 

(6) 



where we have set the Boltzmann constant equal to one. 
From the overdamped Langevin equation , one can de- 
rive a Smoluchowski equation for the probability density 
of the phase a(9, t) (see Ref. [l2| for instance): 



do 
~dt 



1 d 

^Md9 



8U do 
~^9 a + T d9 



(7) 



Notice that the coefficients in this equation are not de- 
pendent on the contact capacitance, as expected for the 
strong damping regime. Floquet's theorem then tells us 
that its general solution has the form 



a(0, t) 



,-At ~ 



(8) 



The exponential factor represents the relaxation towards 
the stationary solution with a time constant A -1 . Note 
that in our case, owing to the time dependence of the 
potential U, the stationary solution oscillates in time 
with the external frequency to. The characteristic fre- 
quency given by lu c — RI c /4>o, where I c = max[I(6)] 
is the critical current, can be considered as the typical 
relaxation rate of the system in the absence of thermal 
fluctuations^. As it is well known for current biased tun- 
nel junctionsii, the observation of well defined Shapiro 
steps in the IV curves requires 10 > lo c . Notice on the 
other hand that for the validity of the adiabatic approxi- 
mation discussed in Sect. ^-B we are requiring ui <§C A^. 
In the rest of this work we will only consider stationary 
solutions (A = 0) and will thus identify o with b. 



B. Current-phase relation 

A crucial ingredient of the model is the current-phase 
relation 1(9). For simplicity we consider a contact with 
one conduction channel of transmission r G [0, 1] (exten- 
sion of the theory to the multichannel case is straightfor- 
ward) . As it is well known from the mesoscopic theory 
of the Josephson effeclii^ the current through the contact 
is carried by the so-called Andreev states with energies 
given by 



e±(9) = ±A SC \ 1 - rsin 



(9) 



(Asc being the superconducting gap) whose separation 
at 9 = it is = 2Asc\/l — t. Note that this Andreev 
gap closes in the ballistic limit r — » 1 (figure [5J). When 
the voltage is much smaller than the Andreev gap (eV <C 
A a), one can assume the system to remain in the state of 
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FIG. 2: Andreev levels as a function of the phase difference 
for r = 0.9 and r = 1. 




decomposition. The step at V = 4>$um/k is then given 
by^ 
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(11) 



where I m denotes the harmonics in the current-phase 
relation, J n are the integer order Bessel functions and 
a = ev ac /hLU. 

When the voltage is comparable to the Andreev gap, 
the adiabatic approximation breaks down since one can 
no longer assume the system to stay in the lowest An- 
dreev level. In the present work we shall restric our anal- 
ysis to the adiabatic approximation and concentrate on 
the effects of phase fluctuations in the current-voltage 
characteristics. 



III. CURRENT- VOLTAGE CHARACTERISTICS 

This section is divided into three parts. In the first 
one we derive the expression of the mean current and 
the mean voltage in terms of the solution of the Smolu- 
chowski equation. In the second one we derive an ex- 
pression for this solution in terms of a matrix continued 
fraction. In the third part we show the numerical results 
for the I — V characteristics. 



FIG. 3: Evolution of the adiabatic current-phase relation 
HUH with the transmission r at zero temperature. 



A. Expressions of the mean current and the mean 
voltage 



lowest energy (this is the adiabatic approximation) and 
obtain the following current-phase relation^: 



1(0) 
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(10) 

When r is small, one recovers the usual sinusoidal rela- 
tion for tunnel junctions with an expression of the critical 
current that matches Ambegaokar-Baratoff formula 1 ^. 
On the other hand, when the contact is strongly trans- 
missive the contribution of higher harmonics can no 
longer be neglected and the current-phase gradually ap- 
proaches the sin (9/2) behavior characteristic of the bal- 
listic limit (figure . These higher harmonics give rise 
to fractional Shapiro steps when a microwave field is ap- 
plied. 

An estimate of the size of these steps at zero tem- 
perature can be obtained assuming that the total bias 
voltage V remains constant. The phase evolution is then 
given by 6(t) — Vt/(po — v ac sin (cot) / ipgLu + 6q. In this 
limit and within the adiabatic approximation one can 
thus obtain an expression for the size of the fractional 
Shapiro steps by introducing this phase evolution into 
the current-phase relation I|1U|) and performing its Fourier 



If we define w = —-^jj [— + T||] , we can rewrite 
@ as: 



da dw 

Tt + 



= o 



(12) 



This equation can be seen as a conservation law for the 
probability (which always holds for stationary solutions) . 
Henceforth, w is a probability current and must be given 
by: 



de v 

w = a— = a— 

at cpo 



(13) 



Let us denote by (...) the mean value of a quantity 
with respect to the phase and < (...) > its mean value 
with respect to time. Our aim is to calculate the current- 
voltage characteristics, that is < 1(6) > as a function of 
< v >, where 

In c oo 



<I(6)> = / d6 dta(9,t)I(8) (14) 

JO J -oo 

1*2-11 poo 

<v> = 0o / d9 dtw(6,t) . (15) 
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FIG. 4: Lattice representation of the Fourier-space Smolu- 
chowski equation. The amplitude of the microwave v ac cou- 
ples the n-chains, while the harmonics of the current-phase 
relation couple the fc-chains. 



J 



The Smoluchowki equation J7J) is periodic in time and 
phase, and so must be the density cr(9, t) and the prob- 
ability current w(0,t). They can thus be expanded in 
double Fourier series: 



a(9,t) = ^ CT ",fc e 



ikO-\-inu)t 



w(9,t) = 22 w n,ke 



ik9-\-inujt 



(16) 
(17) 



The normalization condition for the probability density 
then imposes that 



Cn,0 = #n,o/27T 



(18) 



In Fourier space, the Smoluchowski equation Q reads 



inu>a n k = 



T)M 



-k 2 Ta n ,k - ik(f>oIbO'n,k + ik4>Q ^ Im{<J n ,k-m - VnM+m) 



+k(j)o-^-{a n+ i.k + cr n -i,k) 
ZH 



(19) 
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It is worth noticing that this set of equations can be 
associated with a (non-hermitian) lattice model for par- 
ticles in a square lattice. Each component of the density 
<j n ,k can be associated with a site (n, k) in a Z x Z* 
square lattice. The coupling between chains n and n ± 1 
is proportinal to the ac- voltage v ac , while the coupling 
between chains k and k ± m is proportional to the ru- 
th harmonic of the Josephson current 1(9) (see figure 0J. 
This analogy will be useful for deriving approximate an- 
alytical solutions as discussed in Sect. IIVI 

Note, on the other hand, that temperature appears in 
Eq. Ijl9|l with a factor k 2 . We will see later that in the 
limit of large ac voltage k 2 T actually plays the role of an 
effective temperature for Shapiro steps of order n/k. 

Making use of the orthogonality of circular functions, 
one can easily show that 

< v >= 27r0 o w o ,o = R\Ib - a 0,kl-k) (20) 

fcez 

< W) >= °o,kI-k ■ (21) 



Thus, in order to calculate the current-voltage charac- 
teristics, we only need to know the co.fc's. A Shapiro step 
of order n/k in the IV characteristics is precisely due to 
the presence of a jump in ag^ as a function of the bias 
current. 



B. Recursive solution of the Smoluchowski 
equation 

We now expose how to solve l|19l) numerically in or- 
der to obtain the Fourier components of the probability 
density. 

It is convenient to introduce the vec- 
tors <7„ = (..., a n ,2> Cn.lj er ni _2, •••)) ^ = 
(..., 7_2, J_i, It, I2, •••) and the matrices A„ defined 
by: 
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ikT + 4> Ib Skk 



-Im4>o{&k' ,k-m ~ <5fc',fc+ m ) • (22) 
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We can then rewrite (|19fl in a more compact form, 



we obtain that 



A„S„ = 0O7TFT(<?n+l + OVi-l) + fin,07T~I ■ (23) 



From now on, we take n > 0, and define n = — n. If we 
define the matrices S„ and S„ by 
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FIG. 5: Mean current as a function of the mean voltage across 
the contact, for three different values of the transmission and 
the following parameters (expressed in reduced units as indi- 
cated in the text): R = 10" 4 , T = 10" 2 , v ac = 4.4 x ICT 3 , 
lu/2tt = 12.5 x 10" 3 . 



<i( e j> 



0.14 
0.12 
0.1 
0.08 ' 
0.06 
0.04 ; 
0.02 


-0.02 
-0.04 
-0.06 
-0.08 1 



step 1/3 

step 1/4 

I ,", step 2/5 



step 1/2 



step 2/3 
step 3/5 step 3/4 



FIG. 6: Zoom of FIG. HI 



C. Numerical results 
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A recursive numerical solution of this last equation en- 
ables us to find c?o, and thus the mean voltage < v > 
and the mean current < 1(8) >. The accuracy of the re- 
sults depends on the number of harmonics considered (in 
both, phase and time), that is on the size of the matrices 
(kmax) and the cut-off in the continued fraction (n max ). 
As the temperature is lowered and the ac voltage is in- 
creased the values of n ma x and k rnax required to get a 
good precision increase. The numerical results exposed 
below were obtained with n max = 100 and k max = 50 
which where found to be sufficient to get reliable results 
at the lower temperatures considered. 



In order to illustrate the type of results that are ob- 
tained by the numerical solution of Eq. I|19|) we have 
chosen parameters which roughly correspond to a pos- 
sible experiment based on Al atomic contacts 10 . The 
energy scale in such an experiment is set by the super- 
conducting gap Asc — 180/iteV. Typical microwave fre- 
quencies used in experiments are hut ~ 10~ 2 — 10~ 1 Asc- 
On the other hand, as we stated in the introduction, 
the series resistance should be small compared to the 
resistance quantum Rq = h/Ae 2 in order to neglect 
Coulomb blockade effects at small temperatures 16 . More- 
over, the observability of the fractional steps requires that 
lu > RI c /(f>o, i.e. much larger than the typical relaxation 
rate in Eq. J7J. As I c ~ Asc/fio, this condition implies 
that R < 10~ 2 Rq. Figure shows the current- voltage 
characteristics for different values of the transmission for 
a set of parameters chosen according to this criterium. 
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FIG. 7: Evolution of step 1 with temperature. 



izontal (i.e. labelled by k) chains are nearly decoupled, 
that is when the ac- voltage (resp., the Josephson current) 
is small. Notice that for the typical choice of parameters 
discussed in the previous section RI(9) <C v ac , which cor- 
responds to this last case. The spirit of this approach is 
to decouple the n and k dependences of the Fourier com- 
ponents of the probability density. This decoupling will 
then lead to recurrence relations between the harmon- 
ics that ressemble to those of Bessel and modified Bessel 
functions: 



— Mz) = Jv-i(z) + J„ +1 (z) 
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—X u (z) = l v -x(z) - l v+ i{z) 



(28) 
(29) 
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FIG. 8: Evolution of step 1/2 with temperature. The de- 
crease of the amplitude is clearly much stronger than for step 
1. 



Resistance is measured in units of Rq and all energies 
in units of the superconducting gap A$c (e-<7- m units 
of Asc/e for the voltage and Asc/4>o for the current). 
Naturally, the overall current increases with the trans- 
mission t. As expected, additional Shapiro steps appear 
as t raises. Close to perfect transmission (r = 0.99), one 
can clearly distinguish steps 1, 1/2, 1/3, 2/3, 1/4, 3/4, 
2/5 and 3/5 (figure©. 

The temperature dependence of the I-V characteristic 
is illustrated in Figs. [7\ and [S] showing the behavior of 
steps 1 and 1/2 respectively. As already anticipated, a 
stronger suppression is observed for the 1/2 step. We 
shall analyze the scaling of the fractional steps with tem- 
perature in more detail in the next section. 



Like in Ref. |8|, this analogy can be used to obtain an- 
alytical solutions within this weak interchain coupling 
approximation. 



A. Limit of small ac-voltage and low transmission 

Let us define {Xk',k e Z*} as the solution of the 
Smoluchowski equation l|19|l in the absence of microwaves 
(v ac = u> = 0). In the tunnel limit 1(9) = I c sin#, it obeys 
the equation 

(~ikT + (j)ah)Xk = 0o^(A fc _i - Afc+i) . (30) 

Ivanchenko and Zil'bermanS found its solution by not- 
ing the analogy with the recurrence relation for modified 
Bessel functions T v (z): 



x k = &(-k) 



X k+i%± 



t ) 



4>qIc 

' Q(k) 



(¥0 



To arrive to this expression, one has to impose that the 
probability density is real and normalized (O represents 
the Heavyside function) . If we now introduce a weak cou- 
pling v ac between the n-chains, we can try the following 
ansatz for the Fourier components: a„.k — Kn,fcAfc, with 
K n k — K n k±i- The Kn^'s are then solution of 



In uj4>o 

k V ac 



K n,k — l,k T ^n+l.fc 



(32) 



IV. APPROXIMATE ANALYTICAL 
SOLUTIONS 



Making use of the analogy with the recurrence relation 
of Bessel functions J v (z) , we obtain 



In the quest for analytical expressions for a(9) it is use- 
ful to have in mind the lattice representation of Eq. (|19fl . 
We shall consider two limiting cases corresponding to sit- 
uations in which the vertical (i.e. labelled by n) or hor- 



Kn,k 



(33) 



We thus obtain a generalized expression of the 
Ivanchenko-Zil'berman solution for the normalized prob- 
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FIG. 9: Scaling of the size of Shapiro steps of order n/k with 
temperature for a contact with r = 0.99 at fixed ka — 2.5 
and R = 10 -3 . The full line corresponds to a fit of the form 
1/(1 + ak 2 T). 



ability density in the limit of small a 
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2ir 
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(34) 



(36) 



B. Limit of small Josephson current 

In this limit the coupling between k chains is negligi- 
ble. To zero order in the Josephson current the Fourier 
components of the probability density satisfy the equa- 
tion 

/nuirjM . U r r , T ,\(o) v ac, ( (o) . (o) \ f „ K \ 

Hence is simply given by J n +u k {2ka)/2ir, where 

v k = (-ik^RT + kRI b <po)/(u<f>o)- From and 
we see that, in order to calculate the I — V curves we 

only need a^ k and that it is the only quantity involving 
temperature. Thus, the mean current within this approx- 
imation is given by 



Within this approximation a Shapiro step of order n/k 
arise from the terms in I± k when RI^ ps naj(f>Q/k. The 
above expression suggests that the size of the step should 
scale as J Vk (2ka), i.e. a universal temperature behaviour 
should be observed when the size of the step is plotted 
against k 2 T for fixed ka. 

This approximate scaling is illustrated in Fig. [5] As 
can be observed, the exact numerical results in this limit 
fulfil reasonably this predicted universal behavior. At 
large temperatures this can be well approximated by a 
universal function of the form 1/(1 + ak 2 T), shown as a 
full line in Fig. However, it should be noticed that 
some deviation from universality is always observed due 
to the finite coupling between the fc-chains not included 
in the approximation given by Eq. (|36|) . This deviation 
is more pronounced for the first Shapiro step (k = 1) 
which arises from the larger component of the current 
phase relation. 



V. A MORE REALISTIC CIRCUIT 

Neglecting the capacitance in the usual RCSJ model 
leads to considerable simplification of the circuit equa- 
tions. However this may result in a rather crude descrip- 
tion of the electromagnetic environment in an actual ex- 
periment. We consider in this section a more realistic 
circuit with two resistors and a capacitance as shown in 
Fig. I1UI This model reduces to the previous one when 
r = 0. The presence of two resistors means that there 
are now two sources of thermal noise. We thus have to 
deal with two coupled Langevin equations: 



C 



dV 

~dt 
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h - HO) 
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L R {t) 



(37) 



<po— = V — rl(0) — v ac cosujt — rL r (t) . (38) 
at 

Using standard techniques (see Ref. ^3 for instance), we 
can derive from them a Fokker-Planck equation for the 
probability density W(8, V, t): 



dW 



d_ 



V — v ac cos tot — rl(6) 
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d i RI b - RI{6) - V 



dV 
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T d 2 W rT d 2 W 
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b 2 dv 2 



(39) 
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and noting that in the strong damping equation Q. 
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one easily recovers the Smoluchowski 
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FIG. 10: Generalization of the RCSJ model. 

In order to solve the Fokker-Planck equation J2HJ), we 
need to find on which basis to expand the V dependent 
part of the distribution function W. Following Ref. 0, 
we split the Fokker-Planck operator in a reversible and an 
irreversible part, and note that the T^-dependent term of 
the latter coincides with the Hamiltonian of a harmonic 
oscillator when brought to an hermitian form. This sug- 
gests to consider Hermite polynomials Hk(z). More pre- 
cisely, we will use the following orthonormal basis, 




irMT] 



The distribution will then read 



W(9,V,t)=MV)J2^(V)< 



and the normalization implies 



W nao = 



inujt-\-imQ 



w, 



nmk i 



(41) 



(42) 



The expressions of the mean voltage across the resistance 
R and of the mean current and voltage across the contact 
are then 



« V » 



« 1(9) » = 




-m.O 



(43) 
(44) 



«v» = « V » -r^7 TO W ,-m,o , (45) 



where << (...) >>= J J J dtdVd0(...) represents the 
mean value over the time, the voltage and the phase. 
Thus, again, the knowledge of the n = component of 
the distribution is sufficient to evaluate the quantities of 
interest. 

The irreversible part of the Fokker-Planck operator be- 
ing diagonal in the basis {^k(z);k G N}, we can antic- 
ipate a simple recurrence relation between the compo- 
nents of the distribution in the voltage sub-space, which 
should lead to a continued fraction of matrix continued 
fractions generalizing (|26(l . 

Making use of the orthogonality of circular and ip k 
functions, together with the recurrence relation for Her- 
mite polynomials, 

H k (z) = 2zH k ^(z) - 2{k - l)H k _ 2 (z) , (46) 



(40) we can rewrite the Fokker-Planck equation as 



+imev ac (W„+i, m ,fc + W n -x >m ,k) (47) 



where Exploting its block-tridiagonal structure the set of 

equations (|4*7)l can be evaluated using a recursive algo- 

A k > k i _ _g , I AK n Q e |_ 2t - Tin 2 ) + rithm similar to the one discussed in Sect. III. In order 

to illustrate the effect of this more complex environment 

(1- 

Ak,k— 1 c 
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FIG. 11: Supercurrent peak for a contact with r = 0.9 in 
the circuit of Fig. ITolwith R = 0.001, r = R/W, T = 0.02 
and e 2 /C = 0.0001, 0.0002, 0.0004, 0.0008 and 0.0016 from 
bottom to top. As can be observed the results converge to an 
asympthotic curve corresponding to the RSJ model. 

we concentrate in the analysis of the supercurrent peak 
in the absence of microwaves 18 . Figure ^2 shows the 
current- voltage characteristic around zero bias for a con- 
tact with r = 0.9, r = i?,/10 and different values of the 
capacitance in the circuit. As can be observed the width 
of the supercurrent peak tends to increase as the size of 
the capacitance is reduced. At the same time the height 
of the peak exhibits a slight reduction, which is almost 
negligible at low temperatures. The shape of the super- 
current peak approaches the one found in the standard 
RSJ model when e 2 /C > 2itR. 

VI. CONCLUSION 

We have studied the effect of classical phase diffusion 
on fractional Shapiro steps in quantum point contacts. 



For this purpose we have generalized the standard RCSJ 
model to superconducting contacts with arbitrary trans- 
mission in the presence of a microwave field. In the over- 
damped limit the circuit equations can be mapped into a 
Smoluchowski equation for the probability density of the 
phase difference across the contact. We have presented 
an efficient algorithm for the numerical evaluation of this 
equation. It has been shown that the fractional steps ex- 
hibit a stronger suppression with temperature than the 
interger ones in agreement with preliminary experimen- 
tal findings^. In the limit of large microwave ampli- 
tude Shapiro steps of order n/k exhibit an approximate 
universal behavior as a function of an effective temper- 
ature T e ff = k 2 T and an effective microwave parameter 
Oi e ff — ka. We have also considered the case of a more 
realistic environment including two resistances and a fi- 
nite capacitance for which we derived the corresponding 
Fokker-Plank equation. Our numerical results indicate 
that the main effect of the finite capacitance is to reduce 
the width of the supercurrent peak. We expect that the 
theoretical analysis presented in this work may be useful 
for the proper interpretation of future experiments. 



Acknowledgments 



The authors would like to thank J. Ankerhold, M. 
Chauvin, J.C. Cuevas, D. Esteve, P. Joyez, A. Martin- 
Rodero, P. vom Stein and C. Urbina for interesting dis- 
cussions and suggestions. This work was supported by 
the European Union through the DIENOW Research 
Training Network and by the Spanish CICyT through 
a Picasso contract. 



1 S. Shapiro, Phys. Rev. Lett. 11, 80 (1963). 

2 B. D. Josephson, Phys. Lett 1, 251 (1962). 

3 J. C. Cuevas, J. Heurich, A. Martin-Rodero, A. Levy Yey- 
ati and G. Schon, Phys. Rev. Lett. 88, 157001 (2002). 

4 For a review see, N. Agra'it, A. Levy Yeyati and J.M. van 
Ruitenbeek, Phys. Rep. 377, 81 (2003). 

5 E. Scheer, P. Joyez, D. Esteve, C. Urbina and M.H. 
Devoret, Phys. Rev. Lett., 78, 3535 (1997); E. Scheer, 
N. Agrait, J.C. Cuevas, A. Levy Yeyati, B. Ludoph, A. 
Martin-Rodero, G. Rubio Bollinger, J.M. van Ruitenbeek 
and C. Urbina, Nature 394, 154 (1998). 

6 M. F. Goffman, R. Cron, A. Levy Yeyati, P. Joyez, M. H. 
Devoret, D. Esteve and C. Urbina, Phys. Rev. Lett. 85, 170 
(2000). 

7 M. Tinkham, Introduction to superconductivity (McGraw- 
Hill, New- York, 1996). 



8 Yu. M. Ivanchenko and L. A. Zil'berman, Sov. Phys.-JETP 
28, 1272 (1969). 

9 V. Ambegaokar and B. I. Halperin, Phys. Rev. Lett. 22, 
1364 (1969). 

10 P. Vom Stein, M. Chauvin, R. Duprat, A. Levy Yeyati, 
M. E. Huber, D. Esteve, and C. Urbina, to be published 
in Proceedings of the Vth Rencontres de Moriond in Meso- 
scopic Physics. 

11 P. Russer, J. Appl. Phys. 43, 320 (1972); C. Vanneste et 
al, Phys. Rev. B 31, 4230 (1985). 

12 F. Schwabl, Statistical mechanics, (Springer, Berlin, 2002). 

13 C.W.J. Beenakker in Proceedings of the 14th Taniguchi In- 
ternational Symposium on Transport Phenomena in Meso- 
scopic Systems, edited by H. Fukuyama and T. Ando 
(Springer- Verlag, Berlin, 1992), p. 235. 

14 V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486, 



10 



(1963). 

B. Carvello and C. Urbina, unpublished. 
See for instance G.-L. Ingold and Yu.V. Nazarov in Single 
Charge Tunneling, edited by H. Grabert and M.N. Devoret 
(Plenum Press, New York, 1992). 

H. Risken, The Fokker-Planck equation (Springer- Verlag, 



Berlin 1989). 

For the Shapiro steps we expect a qualitatively similar ef- 
fect due to the finite capacitance as for the supercurrent 
peak. 



